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Abstract 

The growth function of populations is central in biomathematics. The main dogma 
is the existence of density dependence mechanisms, which can be modelled with 
distinct functional forms that depend on the size of the population. One important 
class of regulatory functions is the ^-logistic, which generalises the logistic equation. 
Using this model as a motivation, this paper introduces a simple dynamical refor- 
mulation that generalises many growth functions. The reformulation consists of two 
equations, one for population size, and one for the growth rate. Furthermore, the 
model shows that although population is density-dependent, the dynamics of the 
growth rate does not depend either on population size, nor on the carrying capacity. 
Actually, the growth equation is uncoupled from the population size equation, and 
the model has only two parameters, a Malthusian parameter p and a competition 
coefficient 6. Distinct sign combinations of these parameters reproduce not only the 
family of ^-logistics, but also the van Bertalanffy, Gompertz and Potential Growth 
equations, among other possibilities. It is also shown that, except for two critical 
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points, there is a general size-scaling relation that includes those appearing in the 
most important allometric theories, including the recently proposed Metabolic The- 
ory of Ecology. With this model, several issues of general interest are discussed such 
as the growth of animal population, extinctions, cell growth and allometry, and the 
effect of environment over a population. 
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1 Introduction 



The logistic equati on is a paradigm for populati on biology. T his simple m odel 



Verhulst . 


1838: 


Pearl . 


1927) 



19271 ) or discrete f)Mav 



1976) ver- 



sions describe two fundamental properties of population biology, which are 
(i) the initial exponential rates of growth, and (ii) density-dependent effects, 
like competition under limited resources, indicated by saturation values. The 
discrete logistic equation, in itself opened a new and broad field in biology re- 
l ated to chaot i c behaviours, and for which s ome experimental evidences exist 



( Han ski et al. 



1993 



Gonzalez et al 



2003a ). The continuous version of logis- 



tic growth, although sharing properties with its discrete analog, differs in some 
aspects. It does not show intrinsic bifurcations as the discrete version does, 
and is much more simple to treat analytically. 



Gilpin and Avalal (jl973T ) and iGilpin et al.l (J1976) introduced a model that 
"slightly" generalises the popular logistic equation. Their model, consists on 
modifying the term corresponding to the density- dependence with an exponent 
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9: 



x = px 



X 



where x is population size, the notation x corresponds to time derivative dxj dt, 
and p, 9, and parameters of the model. Compared to the logistic 

equation, their "global model" describes a population that converges in time to 
the same size as the logistic growth, i.e. to the carrying capacity Xoo. However, 
the exponent 9 gives new interpretations to this sigmoid model of growth. If 
9 > 1 then intra-specific competition is high, and the population takes more 
time to reach its asymptotic value, termed carrying capacity. If < 9 < 1 then 
competition is lower and the car rying capacity is reached earlier than in th e 



corresponding logistic dynamics ((Gilpin and Avala 



1973; 



Gilpin et al 



19761 ) 



The ^-logistic model, as it has been termed afterwards, intr oduced a new con 
cept on population ecolo gy that is the ^-selection strategies (jGirpin and Avala 



1973 : 



Gilpin et al. 



1976). Originally, they proposed the model to explain data 



from c ompeting Drosophi la systems after failure to use a Lotka-Volterra-like 



model ( Avala et al 



19731 ). Afterwards, non- competitive versions of the sys- 
tem (i.e. one "allele" or one "species" model) have been used in conserva- 
tion e cology to model avi an population dynamics and calculate extinction 



times ( Saether et al. 



2000), and als o to estimate the effects o f environmental 



stochasticity on population growth (Sa ether and Engen . 



2002J). Other popula- 



tion models have included stochasticity to aid parameter estima tion and study 



the e ffect of environmental changes in caprine populations (jSaether et al 



20021 ). This model h as also been employ e d in 



species abundance (Di serud and Engen 



community ecology to estimate 



2000). The ^-logistic equation is a 



'slightly more complicated model [that] yields significantly more accurate re- 
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suits", using the original words of iGilpin and Avalal (|1973h . 



There are, however, other kinds of regulation terms that have been successfully 
employed to model other kinds of populations and growth. Sigmoid curves in 



particular are attractive 
^-logistic equations. The 



or biologists, but are not necessarily described by 



von Bertalanffvl (|1966l ) equation, for example, is a 



sigmoid curve that is frequently used in allometric modelling, as well as the 
recently pr oposed (and contro versial) curve derived from bioenergetic consid- 



erations by 



West et al 



pertz equation ([Gompertz . 



(200 11). A nother kind of sigmoid is given by the Gom- 



1825I ). which was originally formulated to model 



human demographic data. The G ompertz equat i on ha s become an important 



tool in modelling tumour growth ([Norton et al. 



1976), although applications 



include a wider range of topics. All of these sigmoid share the property of 
reaching carrying capacity, although they have different functional forms (Ta- 
ble 1), which confere distinct dynamical properties: inflection points, critical 
behaviours near x = 0, or rate of convergence to equilibrium. 

However, not all populations obey saturated growth. Among non-saturated 
growths the first classical ex ample is expon ential growth, typically employed 



to describe bacterial cloning 



Hershev 



19391 ) , or simply as descriptors for non- 



regulated conditions of growth. A "general version" of the exponential is po- 
tential growth -which actually shows some kind regulation but do es not reach 



a car rying capacity. Poten tial growth appears in tumour biology ((Hart et al 



9981). early- life evo 



(Calder 



1984: 



in allometry ((Peters 



Rofi, 



986 



ution (Szathmarv and Dcmctcr, 



Dav and 



1983; [Calder, 



Tavl 



1984; 



1997 



1987) 



Ste arnsl. 



Brown and West 



ife hi story theory 



2004), as well as 



20001 ) . Potential 



growth functions is typically a consequence of complex systems where there 
are several levels of organisation having a direct consequence on growth (e.g. 
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fjSzathmarv and Demeterl . 11987c I West et all 119971 )). 



Motivated by the ^-logistic equation, this paper introduces an alternative way 
to interpret and formulate population dynamics models. The description ex- 
plained through out this paper reduces exactly to most common population 
models, including the above-mentioned growth dynamics (resumed in Table 
1). With this new formalism general scaling laws are derived, using initial pop- 
ulation size and carry ing capacities. Th ese scaling laws, include the heuristic 



scaling introduced by 



West et al. 



(200l|) in allometry. 



2 "Mechanics" of Self Regulation 

One of the central issues in population dynamics is to determine the growth 
function that describes a particular population. Growth dynamics in general 
can be expressed in the form 

x = xr(x) . (2) 

The growth rate r(x) is an explicit function of the size of the population x. 
Depending on the nature of the self- regulation, r(x) has different functional 



forms (see table 1.) For a wide review of density depend ence functions, inc. 



ing (mainly) discrete dynamics, the reader can refer to 



Henle et al. 



(200J). 



ud- 



Growth of a population is a contribution of two terms: a replication term, 
whose contribution to the rate is typically a positive constant indicating the 
number of offspring per unit of time, or frequency of duplication, and an- 
other term describing interaction and/or growth inhibition, which is in itself 
dependent on x. 
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Table 1 

Common regulation functions for different population growth models. In all these 
equations x is the size of the population, p is the Malthusian parameter, 6 is the 
competition coefficient, and a is a parameter determined from environmental condi- 
tions. When populations grow to a saturation, a is related to the carrying capacity. 



Model 



Growth 
Rate 



Malthusian Interaction Initial 
Parameter Parameter Rate 



Exponential 



Logistic 



^-Logistic 



Gompertzian 



Potential 



r(x)=a 



r(x)=p(l — ax) 



r(x)=%(l - ax 9 ) 



r(x)=— plog(ax) 



r(x)=ax 6 



p = 



p>0 



p>0 



p>0 



p = 



von Bertalanffy r(x)=—3p(l — ax 1 ^ 3 ) p > 



West et al. (2001) r(x)=-4p(l - ax' 1 ^) p>0 



6 = a (0) = a^0 

6 = 1 a(Q) > 

6 > a(0) > 

6 = a (0) > 

6^0 a(0) ^ 

9 = -i a(0)>0 

9 = -i a(0)>0 



The logistic model, for example, describes growth inhibition with a second 
order term, i.e. x 2 . The ^-logistic generalises this second order term to one of 
an arbitrary order greater than one, expressed by x e+l . In a biological sense, 
this non-linear term is proportional to the frequency of interactions that an in- 
dividual must have in order to produce population growth inhibition (anergy) 
or promotion (synergy). For these examples, the rate will be a decreasing line 
with x for the logistic case, and a decreasing curve of order 9 for the ^-logistic. 



The idea in this paper is to express the dynamics of a population using not 
only its size x as the variable of interest but a decomposition where the rate r 
at which the population grows is considered as a separate state variable. The 
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decomposition requires the knowledge of how r changes in time, thus the start- 



ing point is to determine the dyna mics of r. T 



two d imensional dynamical system f Ermentrout 



lis decomposition results in a 



2002; 



Arrowsmith and Place 



1990) with the variables (x, r) G M + x R describing the population size and 



replication velocity, respectively. 



2.1 Exponential Growth 



In the case of the exponential growth, because r = const = a, the time deriva- 
tive of the rate r = 0. Thus the following trivial system: 



x = xr , 
f = , 



(3a) 
(3b) 



is equivalent to exponential growth. 

Integrating Eq. (3b) gives the constant a, which is determined by the initial 
conditions of the system (x(0),r(0)). Later, substitution into Eq.(3a) recovers 
the original expression in terms of one variable, i.e. x = xa. 

2.2 Logistic Growth 



For the logistic growth, it is necessary to define the new variable r as 



r(x) := p(l — ax) 



(4) 



where p is the Malthusian parameter ((Fisher 



1930T )). and a > is the inverse 



of the carrying capacity. Thus the rate equation again is expressed implicitly 
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as x = xr, and the time derivative for r is: 



r = —pax = —paxr , 

Regrouping, and then summing and subtracting 1 in the parenthesis, it is 
possible to write: 

r = p(l — ax — l)r = (p(l — ax) — p)r . 

The inner parenthesis of the last expression has the explicit form of r. After 
replacing it with Eq. (4), the rate equation becomes: 

f = (r — p)r . (5) 

Therefore, to solve the dynamical system equivalent to the logistic equation 
only one parameter and an initial condition are needed. Actually, the initial 
condition automatically defines the carrying capacity of the population. 

2.3 6 -Logistic Growth. 

The rate for the ^-logistic is defined as 



Following the same methodology as with the logistic growth, it is not difficult 
to demonstrate that the the implicit form for the rate equation is: 



r(x) :=P(l-ax°). 



(6) 



r = 



{Or — p)r . 



(7) 



8 



Although derived from the ^-logistic, this last equation is very general because 
it includes all the growth laws in table 1. In the limit 9 —>■ 1 the logistic equa- 
tion is recovered, and taking jointly the limits 9, p — > Eq. (7) reduces to the 
simple form of exponential growth. Note that from the logistic equation (5) it 
is not possible to take the limit to the exponential formally, since it does not 
show an explicit dependence on 9 (= 1). To recover exponential growth from 
the explicit form of the logistic, the limit would have to be taken as a — > 0. 
But then the rate of the exponential growth will be p instead of a. In the for- 
malism presented in this paper a and p have distinct properties. On the one 
hand, p is defined as a parameter of the system, and as such may have a role in 
bifurcations and global stability, while a is defined fro m the initial conditions , 



so it does not play any role in local or global stability f A rrowsmith and Place 



1990). Also a particular population grows following a predefined replication 
constant p, which is considered to be determined at least in part by intrinsic 
factors, while a is determined extrinsically b y environmental co nditions which 



define the carrying capacity of the system ((MacArthur . 



19621) . Thus in this 



mechanistic interpretation where r determines growth response, the environ- 
mental issues play no dynamical role unless they are explicitly and dynamically 
affecting growth rate. 



3 Generalised Growth Rates 



The three versions of the model studied above, namely exponential, logistic 
and ^-logistic, conform just a part (in fact a minority) of the possible outcomes 
of the system. They were derived by some non negative combinations of pa- 
rameters 9 and p. These, and other dynamics admitting also negative values 
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for 6 and p, conform a dynamical system that generalises most classic types 
of population growth (Table 1). In other words, the model presented herein 
is a unification of several growth dynamics. Resuming, population growth can 
be described in a general form by the two equations: 



referred to from now on as growth equation and rate equation, respectively 

Although the rate equation (8b) indicates that regulation mechanisms are 
independent of population size, per capita response is a contribution of both, 
individual reproduction (related to the parameter p) and interaction with 
other individuals (related to the parameter 6). 

The units of the parameter p are inverse of time (frequency), and it gives the 
characteristic time scale at which individuals down-regulates the reproduction 
rate when, for example, the population approaches an equilibrium state like 
carrying capacity or extinction. The parameter 6, is non-dimensional, but 
it sets the density scale at which the interaction of an individual with the 
population affects its reproduction rate. 

4 Fixed Points and Stability Analysis 

The rate equation (8b) encloses all the information of the fixed points of the 
population dynamics. The rate equation has two fixed points, namely, r = 
and r\ = p/9. Intuitively r corresponds to the non-trivial equilibrium point 



x = xr , 



(8a) 
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of the growth equation, i.e. when the rates become zero the population is in a 
stationary equilibrium between reproduction and mortality. This means that 
r is a steady state under balanced regulation. 

Take for example the explicit form of the rate for the logistic equation 
r = = r(x*) = P - (l - a(x*f) 

that implies x* = a" 1 ^ 9 , and which corresponds to the carrying capacity. In 
this case, the population has a finite size, balanced by reproduction (replication 
at the individual level) and mortality (competition at the population level). 

The biological meaning of the second equilibrium point, r = r 1; is not so 
obvious. From well-known cases, like the logistic equation, it is possible to 
realise that the population has a fixed point in x* = 0. However, the fixed 
point given by the rate equation (6) means that 




which implies directly x* = 0. Thus r* = r± is equivalent to x* = 

The system (8) suggest a third fixed point: (x*,r*) = (0,0). This point, how- 
ever is a paradoxical point, since both, the rate and the population size cannot 
be simultaneously zero, unless p, 9, and a are zero. 



he stability of these fixed points, can be studied with the eigenvalues method 



( Ermentrout . 



2002; 



Arrowsmith and Place . 




. The Jacobian matrix of the 
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system (8) is 



(x, r) 



x 



26r - p 



(9) 



Each of the two eigenvalues of this m atrix, A^ and A r are derived from evaluat 



ing the Jacobian in the fixed points ([Caswell 



20001 ) . For the first fixed point, 



Pq — ( x o, r o) = 0) (where x* is the asymptotic value obtained by equating 
x* = r~ 1 (0), like for example, carrying capacity), leads to the eigenvalues: 



K = 



T 



(10) 



Now, evaluation of the Jacobian matrix in the second fixed point, i.e. Pi 
gives the eigenvalues: 



V = p/e 

A ri = p 



(zi,ri) = (O,p/0) 



A fixed points will be stable if both of its eigenvalues are negative, and unstable 
if any of them eigenvalues is positive. Since p and 9 can take any real value, 
the stability of the fixed points P an d Pi depends only the signs of these two 
parameters, but not in their magnitudes. 

However, some properties are already evident. First, none of the fixed points 
can be foci, since the eigenvalues cannot take imaginary values, thus oscilla- 
tions and cycles cannot occur. Second, the fixed point Pq has a null eigenvalue, 
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\ XQ = (Eq. 10), indicating that there is a invariant set x mv that has "null" 
stability, meaning that the point at which every trajectory intersects x mv , (a) 
depends entirely on the initial conditions (x(0),r(0)), and (b) is fixed. This 
implies that distinct orbits intersect x mv at different points. The dynamical 
equations imply that x mv (also termed stable manifold) corresponds to r = 
for all non-negative values of size x. 

According to the signs of p and 9, Eqs. (10-11) indicate that there are four 
distinct possible sign combinations for the eigenvalues. For each of these com- 
binations, termed regimes, particular patterns in the trajectories occur. This 
suggests several types of equilibria, convergence to equilibrium (which corre- 
sponds to growth dynamics), and transitions between the distinct regimes or 
types of equilibria (bifurcations), as shown in figure 1. 

The first quadrant in figure 1 where p and 9 are both positive, corresponds to 
a phase space that has an unstable node at Pi = p/9. The trajectories that 
start below this point, i.e. r(0) < Pi are attracted to x mv (Fig. 2.) However, 
if the initial conditions are r(0) > Pi, then the orbits are upper unbounded, 
and growth is unlimited. This means that the line r = p/9 is the separatrix 
for the two possible dynamics. 

Maintaining p > and decreasing 9 — > 0, the systems shows a discontinuous 
(i.e. first-order) transition (Fig. 1). Because the position of the fixed point 
Pi = p/9 increases as 9 decreases, at 9 = this Pi disappears at infinity. 
In this case, when 9 = 0, x mv retains its stability (note that the eigenvalues 
associated to x mv , X xo does not depend on 9), and all orbits in the system 
converge to x mv , for all initial conditions. 

When 9 < -and maintaining p > 0- the system is in quadrant II (Fig. 1), 
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Fig. 1. The value of the fixed point r\ (shown as a shade) as a function of the 
parameters p and 9. The four dots inside the four quadrants correspond to the main 
four classes of models, and the four dots on top of the axes are bifurcation points. 
The arrows follow the explanation in section 4. In quadrant I, r\ is an unstable node; 
the point at 9 = 0, following the arrow, is a bifurcation point (n is at infinity). 
In quadrant II, r\ is a saddle. Following the arrow to p = another bifurcation 
point is found. In this point, the invariant set x mv changes its stability from stable 
(attractive) to unstable. In quadrant III, r\ continues to be a saddle. Following the 
arrow to the point at 9 = another bifurcation is found. The stability of t\ changes 
in quadrant IV to a stable node. Finally, returning to quadrant I, another bifurcation 
is found at p = where r\ changes to an unstable node, and x mv becomes stable 
again. The white dashed lines represent the growth rates listed in table 1. 
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Population size x 

Fig. 2. This regime of the phase space corresponds to the first quadrant in figure 1, 
where p, > 0. There are two types of growth. If the initial conditions are below the 
separatrix r(0) < p/6, then the population grows as a ^-logistic. The intersection 
of the orbits at r = correspond to carrying capacities. Initial conditions the are 
above the separatrix (r(0) > p/9), growth in a cooperative way, where interactions 
between individuals accelerate population increase. 

and its properties change: Pi is a saddle (i.e. unstable) point (now there is 
one positive eigenvalue in Eq. 11.) The invariant manifold x mv still retains it 
attracting stability (Fig. 3). 

If the initial conditions r(0) > —p/9, then the orbits intersect x mv as in the 
previous regime (Fig. 3). Note that 9 < 0, thus these initial conditions include 
a region where negative rates diverge (meaning then that population decreases 
in size toward zero). 

From this point in the second quadrant of the parameter space (Fig. 1), de- 
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Population size x 

Fig. 3. In this regime, corresponding to quadrant II in figure 1 with p > and 
6 < two types of growth are possible. The initial conditions above the separatrix 
at r = p/9{< 0) reproduce sigmoid growth curves which converge to a carrying 
capacity that corresponds to the intersection of the orbits at r = 0. If the initial 
conditions are below the separatrix, then the interactions are anergistic and the 
population decreases hyperbolically , and become extinct in a finite time. 

creasing p to zero while maintaining 9 < leads to a continuous (or second 
order) transition. Here the stable varieties x mv and the separatrix Pi collapse 
onto each other. 

Because at this point X Xl = X ri = 0, we cannot infer about the dynamical 
properties of the nullcline at r = with the eigenvalues method. However, 
with a perturbation on each side of the fixed point, it is possible to determine 
the stability. Then if Ar is the perturbation (a trajectory slightly displaced 
from the fixed point), the rate equation for the regime where p = and 6 < 
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is 



Ar = -\6\Ar 2 , 

thus the system always responds diminishing the rate. Consider the solution 
R(t) of the perturbation to the rate equation (note that there is no "first 
order" approximation in the rate equation for this case): 

Ar 

If the perturbation is positive, the rate will be damped to zero asymptotically. 
This means that the population will grow potentially. If the perturbation is 
negative, the rate will decrease to — oo, in a finite time given by t e = (Ar|0|) _1 . 

Thus x mv repels the orbits on its left (initial conditions r(0) < 0, and asymp- 
totically attracts the orbits on its right (initial conditions r(0) > 0). 

In the third quadrant of the parameter space, where 9 < and p < 0, the 
eigenvalue X X1 > but A n < (according to the relationships Eqs. 11). Thus 
Pi is again a saddle point. In this case the trajectories grow x — > oo, and x mv 
is a separatrix that repels the orbits on its neighbourhoods. Figure 4 shows 
several trajectories for this regime. 

Proceeding in direction to the fourth quadrant, maintaining p < and decreas- 
ing 9 — > 0, once again a discontinuous transition is found, where (x,r) — > oo, 
and x mv still repels the trajectories from each side. 

The last regime when 9 > and p < 0, is characterised by being the only 
one where Pi is a stable point (Fig. 5). x rnv is a separatrix: for those initial 
conditions such that r(0) < (i.e. negative rates), the populations decrease to 
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t 1 1 r 




Population size x 

Fig. 4. An invasion-extinction regime, corresponding to quadrant III of figure 1, 
where p, 6 < 0. This regime represents dynamics that show threshold behaviour. 
If the initial rate is r(0) > the population can invade asymptotically-exponential 
with rate p/6. If initial conditions are r(0) < the population may go extinct in 
finite time with the rate decreasing hyperbolically. 

zero sigmoidally, and if the initial conditions are r(0) > 0, then the population 
increases hyperbolically to x — > oo. 

5 Carrying Capacity 

In the growth dynamics, the explicit rate equations (Table 1) involve a con- 
stant a which does not appear in the implicit form of the dynamical system 
(8). Integrating the growth equation leads naturally to the constant a, as it 
did in the calculations of the exponential and ^-logistic equations in section 
2. Once the initial conditions are denned the constant a is determined. In 
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Population size x 

Fig. 5. An invasion-extinction regime corresponding to quadrant IV in figure 1 with 
p < 0, 8 > 0. This regime is of threshold type. Invasion occurs in a synergistic way 
when r(0) > 0, although when the conditions are met for extinction, i.e. r(0) < 
the population vanishes decreasing exponentially. 



quadrants I and II x mv is stable, and the point at which the orbits intersect 
x %nv correS p 0n d to the carrying capacity of the system, x^. The value of a 
is related to the carrying capacity ) I . Moreover, all sigmoid 

dynamics are represented in these two quadrants. When 9 > the trajecto- 
ries correspond to logist ic and ^-logistic dynami cs, whose app lications range 

1973) to mammals (paether et al 



from populations of flies f G ilpin and Avala 
2002). For 6 = the dynamics correspond to Gompertzian growth, which has 
been widely applied in t umour biology to investigate distinct aspects of tumour 



response and regression (Norton et al 



as microbiological models 



1976 



Nor 



Kozusko and Baizer 



;on and Simon 



19771 ). as well 



2003). Whenever 9 < onto- 



genetic growth laws, like the model bv Ivon Bertalanffvl (|1957f ). or the model 
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bv lWest et al.l (|2001l ). are included in this regime. 



Thus by itself a plays no particular role in the stability of the system. In this 
formalism the growth rate is an adaptive mechanism that responds to the en- 
vironmental conditions, which enter as the initial conditions. Different initial 
conditions result in different growth models. Previous formulations included 
the carrying capacity as a predefined constant, entering explicitly in the model 
as the parameter a or x^. This constrains the dynamics to a particular model, 
and also s uggests that the g rowth of a population is not dependent of the envi- 



ronment ((Mac Arthur 



1962). In this formalism the growth rate is an adaptive 



mechanism that responds to the environmental conditions. The same model 
can grow toward distinct carrying capacities (or even change its qualitative 
form to an unlimited growth) just by changing the initial conditions. Explicit 
models with predefined carrying capacity, need extra assumptions for model 
selection. In this formulation the mechanism for regulation remains robust 
against environmental changes. 

6 Scaling Laws 



Several systems that look very different may actually have exactly the same dy- 
namics when they ar e properly scaled in time and size. Scaling is common phe- 

" x]). However, in biology only recently scal- 



nomenon in physics ((Kadanofi, 
ing hypotheses in dynamical processes are being proposed flWest et al 



Brown and West 



2000; 



West et al. 



2001 



Brown et al. 



2004; 



1997 



Kozlowski and Konarzewski 



2004). Whenever systems can be scaled, it is no longer important (in the dy- 
namical sense) the precise values of the constants. The growth law of the 
population remains invariant whenever proper relationships between the con- 
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stants and the variables are considered. 

To determine the scaling laws for the solutions of the system (8) it is convenient 
to replace r — > Or, which gives 



x = -rx , (12a) 

o 

r = r{r — p) . (12b) 



Denoting the solution for the system (8) as (X(t), R(t)) growth equation yields 
X(t) := x(0) exp ^ J R(s)ds^j , (13) 

where x(0) is the initial condition for x(t). Rearranging this system, we get 

fx(t)\ 



x(0) 



exp J R{s)ds\ . (14) 



The right-hand side of the last equation is independent of 6. Thus populations 
described by the system (8) are always scalable to their initial sizes, and 
interaction exponent 9 (except, maybe for the bifurcation points.) 

Further rescaling is possible for the right-hand side. The solution to the rate 

equation can be written in the form 

-i 



R(t) := p 



1 + e 



P t,_P_ 



(15) 



r(0) J 
and changing the time variable as 



^- ]ot W)-T (16) 



and also changing properly the differential in the integral in Eq. (13) to dT' 
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pdt, then the result of the integral is, in scales of T 



'xyy 



l + e 
1 + e 5 



r 



(17) 



The variable T depends on the initial condition r(0) whose meaning is not so 
obvious. 



Although scaling laws are general, in practice (e.g. in allometry) scaling laws 
are relevant when populations reach a carrying capacity. Suppose then that 
the population achieves a carrying capacity x^, given by x^ = cT 1 ! 6 '. Using 
the explicit form of r(t) (Eq. 6), it results that the term inside the logarithm 
in the rescaled time (Eq. 16) is: 



P 



'x(oy 



r(0) 



1 . 



(18) 



Substituting the relation (18) into the time transformation (16), and rearrang- 
ing terms in equation (17): 



'x(t)' 



1 - 



'x(oy 



- 1 



(19) 



To conclude, define the new scaled variables as: 



T = pt — log 



\ J-oo / 



(20) 



with which the general scaling law obeys 



X = 1 - e" 



(21) 



This scaled dynamical expression is valid for all (non zero) sign combinations 
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of p and 6. 

6.1 Scaling at the bifurcation points 

The scaling law (Eq. s21) cannot be defined for Gompertzian or potential 
growth, because the general solutions for the size and rate equations (8) are 
not necessarily equations (13) and (15). It is necessary to proceed calculating 
the solutions of Gompertzian and potential growth, and then derive the scaling 
laws. 

6.1.1 Potential Growth 

The solution for the potential growth dynamics is given also by Eq. (13). 
Consider then, that the solution to its rate equation is 



In this case, the initial condition r(0) is not expressed in terms of a carry- 
ing capacity. Potential growth or decay does not reach an asymptotic value. 
Thus r(0) can only be defined in terms of the integration constant a. Upon 
integration and rearranging of terms, the following form is found: 



R(t) :=r(0) (l-r(O)t) 



i 



(22) 




(23) 



The scaled variables can then be defined as 




(24) 



< 



t = at 
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with which the general scaling law obeys 



X = 1 -r , 



(25) 



that is simply a decreasing line. 



6.1.2 Gompertzian Growth 



Consider now the solution to the size equation for the Gompertzian Growth 



x(0) 



= exp J R{t')dt' . 



(26) 



The main difference between Eq. (13) and the last equation is that the former 
can be scaled with the exponent 9, which for the Gompertzian Growth case 
is zero. 



Consider then the solution R(t) to the rate equation: 



R(t) := r(0)e- pt 



(27) 



Thus the general solution for Gompertzian Growth is 



x(0) 



exp 



log 



'ar(O)' 



(e-* - 1) 



(28) 



where the initial condition r(0)is expressed as: 
x(0) 



r(0) = -plog 



(29) 



Eq. (28) can be rearranged to give 



x(t) 
log = exp 



X r 



-pt + log log 



X c 



x(0) 



(30) 
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Thus, denning the dimensionless variables: 

X = log(Xoo/x(t)) 

(31) 

r = pt- loglog(xoo/x(0)) 

then the scaled dynamics results as 

X = e- T , (32) 

that is a decreasing exponential. 

Although the scaling behaviour is completely defined for all the dynamics in 
the phase space, it is still interesting to determine how the scaling law Eq. 
(21) is related to the scaling laws Eqs. (25,32). 

Direct evaluation of p = or 9 = does not give the scaling laws for Gom- 
pertzian or potential growth. However, Taylor expansion on the parameters p 
and 9 to the linear term around the critical values, result in the scaling laws 
of the critical points. 

Note that in order to give a precise meaning to the transformations Eqs. (20), 
the initial condition for the rate equation was transformed to its explicit form. 
It is necessary, however, to make this transformation after the limits are taken 
in the scaled variables Eqs. s(20). 

7 Discussion 

The simple model derived in this paper is rich in qualitative solutions since 
it resumes several growth rates that often appear in the literature, which in- 
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elude several levels of biological organisation. The examples alluded in the 
text range from cellular populations of procariots, cellular populations in eu- 
cariots, in ontogeny and cancer, to population biology of mammals and birds, 
to community ecology. It is a nice result that all of these kinds of growth can 
be described by such simple equations that resumes the main features of pop- 
ulations, in the traditional sense of showing density dependence, and in the 
distinct interpretations introduced in this paper. 



7. 1 Extinctions and Invasions 



Environmental changes, are know to "unbalance" some populations. In this 
model, environment enters as initial conditions; translations of the r-component 
in the phase space can model changes in the environment. Suppose that a pop- 
ulation is in its carrying capacity, and suddenly the trajectory is perturbed. 
This will placed the trajectory in an orbit out of equilibrium. Thus the popu- 
lation will converge to a new carrying capacity. However, if the perturbation 
is strong enough, then the orbit where the dynamics is placed could belong to 
a basin of attraction leading to extinction or invasion. 

As a first ex ample, consider glo bal warming. The metabolic theory of ecology, 



proposed by 



Brown et al. 



( 20041 ) proposed considers that the carrying capacity 
of a population is temperature-dependent through a Boltzmann factor. This 
dependence can be expressed as x^ = K exp(E/kT), where K is a parameter 
depending on mass, resources, etc., E is the energy of the limiting metabolic 
reaction, k is Boltzmann's constant, and T the absolute temperature. 

Suppose a population that is in its carrying capacity at a temperature Tq 
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Zoo = Koexp(E/kTo). If suddenly the temperature increases to Ti, the change 
in the rate due to temperature increase can be calculated, introducing the 
Boltzmann factors for population size in the explicit form of rate (e.g. Eq. 6): 



Ar 



P 




exp 



' 9E 



< 



(33) 



This means that the population is taken out from equilibrium. Carrying ca- 
pacities exist whenever p > (for any value of 6). According to Eq. (33) 
T he population will a ttain a new lower carrying capacity. Thus the model 



of 



Brown et al 



(20041) and this formulation imply that temperature cannot 
induce population extinction. To induce an extinction the exponential term 
would have to change sign, which is not possible for any temperature (it is 
always possible however, to consider such changes that although they theo- 
retically do not imply extinctions, numerically are so small that in real life 
populations could disappear.) 

Kin selection, is a second example. The saturated dynamics, comprised in 
quadrants I and II of figure 1, have other co-existing behaviours which support 
Hamilton's rule . which points out that if the cost of an altruistic behaviour is 
such tha t it benefits a genetica 



selected ([Hamilton . 



1963 



ly related individuals, then the strategy can be 
1964J). In this way, a population consisting of coop- 
erative individuals can grow faster than expected by exponential models. This 
can happen when a < in quadrant I (Fig. 2), because growth is accelerated 
as a result of the interaction between the individuals. Notice that for a' < p/6 



a\x > a 



(34) 



for all x] and even if a' > p/9, there exists some x such that for x > x the 
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inequality holds. Thus in this region of the phase space, population growth is 
faster than exponential. 



However, Hamilton's rule has another consequence. This is that aggressive 
behaviours can also be selected, provided that damage is induced to "nega- 
tively related" (i.e. unre 



Gardner and West 



ated) individuals (as it has been reported for wasps 
( 20041 ).) This kind of behaviour is recovered for quadrant II 
(Fig. 3) where any interaction between individuals results in mutual annihila- 
tion indicated by a > 0s. When the initial conditions are r(0) < —p/9, the rate 
decreases. Solving the rate equation (8b), shows that r — > — oo asymptotically 
(Fig. 3) when t — *> t e , where 



log 



r(x ) 



r(x ) 



thus extinction occurs in finite time. 



Another example, m entioned in the text above, is the extinction of sparrows 



(Saether et al 



2000J). In order to have the risk of finite-time extinctions from 
previously stable populations, it is necessary that (a) the population dynamics 
belongs to quadrant II of figure 1, and (b) there is a perturbation such as 
mentioned above. The estimated mean value for the population of sparrows 
is 9 ~ 1, indicating that the population is logistic. However, the estimate 
distribution for 9 allows a small but not negligible probability for —1.5 < 9 < 
0. If this is the case, then a real risk of finite time extinction exists. 
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7.2 Life histories 



In life history theory, the central problem is resource allocation for adap- 
tive strategies. Survival and reproduction in distinct stages determ ine of the 



net growth of an individual (|Dav and Tavlor . 



1997 



Stearns 



2004J). In non- 



reproductive stages of life, i.e. before maturity, resources are mainly devoted 
to growth. The energetic content allocated to growth depends on the size of 
the body. This dependence follows a potential growth function, where the 
exponent 6 is indicative of s ome length scale of the physiological processes 



( Calder 



1984; 



Stearns 



2004). 



According to the rate equation (8b), initial stages in development have to 
be dominated by the term 8r 2 , which implies that resources are devoted to 



growth. This term is dominan t at low densities, whe n th e exponent 6 [< C 



n r 



West et al 



(200 



mi) 



) 



This condition is met in in the lvon Bertalanffvl ([19571 ) and 
equations, and well as in most allometric growth relationships ((Calder 
When reproduction becomes a priority, less energy is devoted to growth. In this 
case the term pr is not negligible, indicating that the organism is partitioning 
the resources between growth and survival. 



7.3 Allometry and scaling 



The scaling laws derived in section 6 are general formulas showing that scal- 
ing is rather a rule than an ex ception. Th i s give s an broader view to the 



open discussion of whether the 



West et al. 



not (Kozlowski and Konarzewski 



( 20011 1 equation is legitimate or 



20041 ). In terms of these formulations, al- 



though numerically different, West's equations and its classical competitor, 



29 



van Bertalanffy's equation, have the same qualitative behaviour. 



The model (8) reproduces the lvon Bertalanffvl (|l966T ) equation. This equation 
can be written in the following form: 



x = ax 



(35) 



It is possible to rearrange this equation, to express it as the system (8), for 
which the parameters are then: 



P 



3 " 



1 

3 ' 



(36) 



The exponent —1/3 follows from the hypothesis that mass is proportional to 
the third power of length, and the parameter b is related to individual repro- 
duction. The parameter a is related to the carrying capacity of the population, 
thus it does not appear in the transformations. 



Anot her example in this regime is West's ontogenetic growth equation ([West et al 
20011 ). given by 



x = am 



1/4 



1 - 



x 



1/4' 



(37) 



The parameters for the rate equation are then: 



p = l^ 1 ' 6 

L.I, -v-. 



1 

4 ' 



(38) 



The exponent —1/4 derived from fractal patt erns of fluid transport systems 



like c irculatory system or plant vascularisation dWest et al 



20041 ) and is also supported by empirical data ([West et al 



997; 



20011 ) 



Brown et al. 



The parameter p in this case depends explicitly on the carrying capacity x c 
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This suggests that in ontogenetic growth, these two quantities could be corre- 
lated, implying that the macroscopic growth i.e. the cell population, is "trans- 
mitting" information to the microlevel, i.e. single-cell dynamics and thus it 
could imply existenc e of self-organisation (this has been demonstrated for 



Gompertzian growth ((Molski and Konarski . 



2003). 



In both models, van Bertalanffy and West, p > and 9 < 0, belonging to the 
same region in the parameter space (Fig. 1) 



7.4 Cancer 



The growth of tumours have been very w ell studied throu gh distinct kinds of 



mathematical and computational models ((Wheldon 



1988). In vitro and exper- 



iment al frameworks have shown that they grow according to a Gompertzian 



Law ([Norton et al. 



19761 ). The Gompertzian law, a sigmoid curve that grows 
toward a carrying capacity, implies that resources in a tumour are devoted to 
reproduction, as indicated by 9 = 0. In the rate-based scheme the orbits are 
attracted to the invariant set x mv for any initial condition. 



Molski and Konarskil (j2003l ) demonstrated that the Gompertz equation is the 
result of self-organisation (cooperativity), in such a way that the response 
of each cell is coherently correlated to t he state of the whole population. In 



their model. 



Molski and Konarskil (J2003J) found a Gompertzian regression rate 



(shrinking) solution. In the formalism presented inhere the Gompertzian re- 
gression solution corresponds to p < and 9 < with initial conditions r(0) < 
0. This kind of decreasing size tumours appear under external perturbations 
(which actually leave the rate equation invariant) are applied, like radiation or 
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chem otherapy treatments ([Gonzalez et all I2003at IP. de Vladar and Gonzalez , 



8 Concluding Remarks 

The ^-logistic equation has become a paradigm in ecology. Modelling popula- 
tions with it has been an important tool to confront actual problems about 
density-dependent ecology. The transformation introduced in this paper gives 
a good insight into the meaning of the quantities appearing in the equations, 
namely p, 8, and a (either in its interpretations as carrying capacity or not). 

There are, of course, more examples for each of the growth types described 
in this paper. However, more interesting is that there are regimes that have 
not been reported. This is not surprising, because they conform distinct types 
of indeterminate growth, which usually is assumed to be "exponential" . How- 
ever, these distinct types of explosions, can have important consequences in 
disciplines like biotechnology, where a strict control of growth is necessary. If 
for some reason, a population is wrongly manipulated, such that it spreads 
"indeterminately" then the distinct types of growths should be managed dis- 
tinctly. 

However, in the opinion of the author, the most important result is that the 
rate equation is explicitly independent of the population size. The results pre- 
sented in this paper, are derived from a simple mathematical transformation, 
which surprisingly results in a very broad class of regulatory mechanism. Al- 
though this is a result that may apply only to the simple systems included in 
this work, it is puzzling why and how the regulatory mechanisms act. 
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The two terms of the rate equation (8b), from a more abstract perspective, 
correspond to two processes that constitute regulation: reproduction, which 
comes from an individual level, and sensitivity to population interaction. De- 
pending on the context, the sensitivity to the population can be of synergistic 
or anergistic nature. The two regulation processes, could be thought of as 
fragmentation and condensation reactions: 

x — 2x (reproduction) 

(39) 

Nx — > Mx (population sensitivity) 

If M > N then the population experiences synergy in growth (i.e. population 
interaction promotes growth). If M < N then the population experiences 
anergy in growth (i.e. population interactions avoid growth). 

The term —pr in the rate equation indicates that the population is growing 
or "relaxing" to a fixed point. Since p is the inverse of the relaxation time for 
the rate, then the bigger p is the smaller the time to let the mechanism to 
relax, and thus the fastest to reach limiting population size at x mv . 

The second term of the rate equation is related to the interaction between 
individuals in the population. This term can be compared to a "potential" 
indicating some kind of resource potentiating (either synergistically or aner- 
gistically) from the interaction. The kind of interaction, is given by the sign 
of 9 and by the environmental conditions, i.e. by a. 

The relationship between p and 9 determining the distinct types of growth 
rates, gives distinct types of behaviours for distinct initial conditions which can 
result in cooperative, competitive, or aggressive strategies. These are strate- 
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gies that can be sought directly from the rate equation. If a population is 
behaving cooperatively, then it means that the rate begins over a threshold 
such that it grows unlimited, because there is a benefit improving growth re- 
sulting from the interaction. In this case, the resources have to be unlimited, 
so cooperativity improves resource allocation for reproduction. In the case of 
competition, usually the scenario is that resources are limited, and there must 
be an equilibrium between reproduction and survival. But if the population 
presents aggressive behaviours, then the initial rate of the population is below 
the threshold where it goes extinct in finite time or exponentially. 

The exponential growth is a particular case in which no regulation mechanism 
(thus interactions) is present in the population. Increase is based only in indi- 
vidual reproduction, and the per capita response is totally independent of the 
state of the system. 

The distinct qualitative solutions (regimes) for population growth have an 
underlying symmetry. The regression equation for population shrinking is in 
general obtained by inverting the time arrow, changing t — > —if. However, 
time-inverting the rate equation does not produce the desired result. In order 
to obtain the regression dynamics from the size-rate decomposition, besides 
inverting the time arrow, it is necessary to invert the rate variable r — > —r'. 
Thus a time reversed equation results in the transformed system 

x = xr' , (40a) 




Notice that these transformations are equivalent to change the sign of the 
Malthusian parameter p. This transformation for population shrinking, can 




{Or' + p)r . 



(40b) 
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be regarded as a reflexion of the parameter space (Fig. 1) with respect to 
the (9-axis. Thus, although most of the relevant population dynamics are in 
quadrants I and II in figure 1, their corresponding regression dynamics are in 
quadrants III and IV. 

The decomposition presented in this work, is a change in the paradigm of 
population dynamics. The Eqs. (8) are very general, but still simplistic because 
there are ma ny biological aspects left aside. Take for example the Allee effect 



(jAllei 



1931). This is a density dependent growth mechanism, not represented 
in the rate equation. Actually, including Allee effect in population growth, 
leads to a polynomial rate equation of higher order than 2. There is however 
no general law that can be derived (at least, including the forms found in 
the literature). Other example in non-random mating, which could lead to 
distinct forms of density dependence, maybe through a frequency-dependent 
parameter 9. In this case the form of the rate equation derived here will not 
necessary be valid. 

It remains to investigate based on life history theory for which kind of resource 
allocation the density-independent rate equation can be derived. This is a work 
currently under development that is expected to help to give conclusions about 
other biologically relevant aspects not included in this work. 
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